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Limited Imitation Contagion on Random Networks: Chaos, Universality, and 
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We study a family of binary state, socially-inspired contagion models which incorporate imitation 
limited by an aversion to complete conformity. We uncover rich behavior in our models whether 
operating with probabilistic or deterministic individual response functions, both on dynamic or 
fixed random networks. In particular, we find significant variation in the limiting behavior of a 
population's infected fraction, ranging from steady-state to chaotic. We show that period doubling 
arises as we increase the average node degree, and that the universality class of this well known 
route to chaos depends on the interaction structure of random networks rather than the microscopic 
behavior of individual nodes. We find that increasing the fixedness of the system tends to stabilize 
the infected fraction, yet disjoint, multiple equilibria are possible depending solely on the choice of 
the initially infected node. 

PACS numbers: 89.65.-s,87.19.Xx,87.23.Ge,05.45.-a 
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The structure and dynamics of real, complex networks 
remains an open area of great research interest, partic- 
ularly in the realm of evolutionary processes acting on 
and within networked systems 16]. Here, motivated 
by considerations of social contagion — the spreading of 
ideas and behaviors between people through social net- 
works and media — we explore an idealized, binary-state 
social contagion model in which individuals choose to be 
like others but only up to a point: they do not want to be 
like everyone else [7Hl0j. We term such behavior 'Lim- 
ited Imitation Contagion.' We build naturally on pre- 
vious studies of threshold models of contagion [Ill4l4| . 
and our model can also be seen as a specific subfamily 
of dynamical Boolean network models [3, [l6| . We show 
how macroscopic network structure overrides microscopic 
details, and we find complex dynamics whose character 
moves from universal and predictable to particular and 
unpredictable as we allow the system to become increas- 
ingly deterministic. 

In constructing our model, our main interest is in 
understanding how spreading by Limited Imitation Con- 
tagion on random networks behaves under three main 
tunable conditions: (1) Social awareness: the rate of 
contact between individuals; (2) Social variability: the 
extent to which friendships are fixed; and (3) Social 
influence: the character of individuals' responses to the 
behavior of others. 



To begin with, wc consider a binary state model 17 1 
for which individuals are either in a base state Sq or an 
alternate state S±. We assume individuals interact over 
an uncorrelated random network (we treat correlated net- 
works in 18(), which may be dynamic or fixed. For sim- 
plicity, and due to the richness of the dynamics we find, 
we employ standard Erdos-Renyi networks which possess 
Poisson degree distributions. We take time to be discrete 
(t = 0, 1, 2, . . .), and we prescribe each node's degree k 



at t = 0. In a dynamic network, when node i updates, 
it samples the states of ki randomly chosen nodes (i.e., 
the system is a random mixing model with non-uniform 
contact rates). For a fixed network, node i repeated- 
ly samples the same ki nodes. We further restrict our 
attention to single-seed contagion processes wherein all 
nodes are in state 5*0 at time t = 0, with one randomly 
chosen node in state S\. 

The contagion process is manifested through the 
response functions of individual nodes. We allow nodes to 
update synchronously (in 18|, we consider a fourth con- 
dition of social synchrony — the degree of timing coher- 
ence between individual decisions). Node i's response 
function Fi : [0, 1] M> [0, 1] gives the probability that node 
i will be in state Si upon updating, where the argument 
taken by Fi is the fraction of nodes sampled by node i 
that are currently in state Si, (f>ij. We investigate two 
kinds of response functions, probabilistic and determin- 
istic, both of which incorporate the characteristic of the 
adoption probability growing with low and then dimin- 
ishing with high values of perceived popularity. 

In Fig. [T]4., we show an example of a probabilistic 
response function, the tent map, which is defined as 
T r (x) = rx for < x < | and r(l — x) for | < x < 1. 
Taking the r = 2 case — for which the iterative map of 
the interval is fully chaotic and the invariant density x 
is uniform on [0, 1] — node i adopts state <Si with proba- 
bility T2(<fiij)- We use the tent map T 2 for several rea- 
sons: (1) As a standard iterative map of the interval, 
the tent map's dynamics are both interesting and well 
understood [19|; (2) The tent map captures a probabilis- 
tic flavor of the adopt-when-novel, drop- when-ubiquitous 
behavior we aim to model; and (3) We can construct a 
simple and elegant connection with deterministic conta- 
gion processes which we describe next. 

In Fig. [Tj3, we show an example of a deterministic 
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FIG. 1: Examples of probabilistic and deterministic response 
functions capturing Limited Imitation Contagion dynamics. 
At each time t, nodes use their given response functions to 
update their own state based on the perceived fraction of their 
neighbors in state Si, cj)i,t- We construct the tent map T2 
(see main text) shown in (A) by averaging over deterministic 
response functions of the kind shown in (B) by considering a 
family of the latter with 'on' and 'off' thresholds uniformly 
distributed in [0, |] and [3,!] respectively. We then build 
networked systems whose macroscopic character is tent map- 
like but differ strongly at the microscopic level. 

response function which is characterized by 'on' and 'off' 
thresholds, cf> on and o g. Node i will only adopt or 
remain in state Si if it perceives the fraction of oth- 
ers in state Si to lie between its on and off thresholds: 

4>i,on < 4>i,t < <f>i,o&- 

Now, for these deterministic 'on-off' response func- 
tions, we examine the special case where <p on is distribut- 
ed uniformly on [0, ^] and <p s likewise on [i, 1]. Aver- 
aging over all deterministic response functions created 
in this way, we obtain precisely the tent map T%. Col- 
lectively, we then have the same on-average behavior in 
both the probabilistic and deterministic cases, and this 
allows us to profitably explore the effect of varying the 
specific response dynamics at the micro level as well as 
interaction patterns. 

We now examine the behavior of four subfamilies of our 
model that vary in terms of node response functions being 
probabilistic or deterministic (P or D), and whether or 
not network connections randomly rewire or are fixed (R 
or F). We refer to these model classes as P-R, P-F, D- 
R, and D-F. Each class is indexed by the parameter of 
average node degree k avg , and, as per our design, all four 
systems have the same on-average response functions (the 
tent map T2) and degree distribution (Poisson). 

Our first focus is on the long term behavior of these 
four networked systems. Key to our understanding are 
the well-known behaviors of the tent and logistic maps 
acting iteratively on the unit interval, the former given 
above and the latter by L r (x) = rx(l — x). We reproduce 
the orbit diagrams for these models in Figs. [2]A_ and B. 
Both systems are controlled by amplitude parameter r, 
whose increase leads to changes in their invariant den- 
sities, famously resulting in the bifurcation diagram for 



the logistic map. The two maps, while topologically con- 
jugate, produce distinct bifurcation diagrams. 

As we show in Fig. Hp-F, the signature of the micro- 
scopic response function of the tent map is erased by the 
network dynamics of the four model classes. Macroscopi- 
cally, we see four orbit diagrams analogous to the logistic 
map's characteristic bifurcation diagram, corresponding 
to Fig. [2j3 rather than Fig. [2]A_. We see that increasing 
the average connectivity of the network fc avg is equivalent 
to increasing the logistic map's amplitude parameter r, 
and the system moves along the period-doubling route to 
chaos. 

While appearing to belong to the same universality 
class, the orbit diagrams of the four models differ impor- 
tantly in detail, most profoundly for the fully determin- 
istic D-F class. First, we observe the four systems pro- 
duce three distinguishable orbit diagrams, since for large 
enough systems, the two random mixing classes P-R and 
D-R must exhibit the same macroscopic behavior. For 
the D-R class, random rewiring overwhelms the fact of 
each node having a fixed deterministic response function. 

We are able to comprehensively explore the P-R and 
D-R models via both simulation and analysis, finding 
excellent agreement. We show only simulation results 
here, making full comparisons in [lij . and record the core 
analytic description below by directly following the P-R 
model. The infected fraction 4> t (equivalently, the prob- 
ability that a random node is infected at time t) is given 
by 4>t+i = J2kLo PkGk(pt), where P k is the probability a 
node samples the states of k randomly chosen nodes; Gk 
is the probability that such a node is active at time step 
t; and pt is the probability that a randomly chosen edge 
emanates from a node infected at time t. The dynam- 
ics are encoded in the update map for the edge infection 
probability: 



Pt+i = 




where G k {p t ) = £ } =0 - p t ) k ^f (j/k) , and 

where f(x) is the node response function (here, the tent 
map) . In [18| , we show that in the limit of infinite average 
degree /c avg and providing the width of the degree distri- 
bution grows sublincarly with respect to fc avg , the dynam- 
ics approach that of the response function /. However, in 
the approach to that limit, the form of the {Gk] induces 
a quadratic maximum for the edge infection update map, 
Eq. ([T]), and the system thereby shares the universality 
class of the logistic map. Thus, in terms of the emer- 
gent, global behavior, we see that the interaction pat- 
tern of the network overrides and masks the microscop- 
ic response function. More generally, we find the same 
universal behavior for any unimodal, concave response 
function flil ]. 

Moving to the fixed network models P-F and D-F, we 
find that these are less amenable to direct analysis, and 
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FIG. 2: Comparison of bifurcation diagrams for networked systems with Limited Imitation Contagion based on tent-map 
response functions. (A) and (B): For reference, standard bifurcation diagrams for the tent map and the logistic map operating 
as maps of the unit interval, showing distinct universality classes. The vertical axis is normalized invariant density \- (C-F): 
Orbit diagrams for Limited Imitation Contagion acting on standard random networks for four model classes with probabilistic 
or deterministic responses and network interactions as indicated. The tunable parameter is average degree fcavg- In F, the 
vertical line indicates fc avg = 30, the system examined in Fig. [3] Simulation details: network size JV=10 4 ; one random seed per 
network; granularity of 0.1 in fc avg ; binning size of 0.001 in (f>; for P-R, D-R, and P-F systems, we ran simulations for 10 4 time 
steps on 10 3 networks, ignoring the first 10 3 values of cj>t\ and for D-F systems, which exhibit long transient behavior (see text), 
we ran simulations for 10 5 time steps on 200 networks, ignoring the first 5 x 10 4 time steps. 



report our findings for simulation results only. Consider- 
ing the P-F model first, we see in Fig. [2fi! that the orbit 
diagram has slightly moved to the right — the bifurcation 
points now occur for slightly higher values of fc avg . Thus, 
the fixedness of the network appears to induce some mod- 
est changes in the dynamics, though we cannot discount 
finite size effects. 

Our final model class, D-F — to whose description we 
devote the remainder of the paper — is the most struc- 
tured, and it generates considerably distinct and intricate 
behavior. Once the network and set of response functions 
is realized, the system is now completely deterministic, 
and we find surprising changes in the system's behavior 
at both macroscopic and microscopic scales. In Fig. [2f\ 
we see that the orbit diagram for the D-F model col- 
lapses abruptly after several rounds of period doubling, 
which are themselves relatively compressed. Our simula- 
tions suggest that above an average degree of fc avg — 16, 
the macroscopic dynamics always collapse. The collapse 
appears to favor a fixed point macroscopic state, around 
<p = 2/3, which is the fixed point of the tent map Ti. 
However, a closer examination of the D-F class's poten- 
tial dynamics reveals a far more subtle story. 

In Fig. |31 we summarize the possible dynamics of a lone 
realized network (7V=10 4 , fc avg =30) with a single set of 



fixed deterministic response functions (we explore other 
values of fc avg in We exhaustively ran N = 10 4 

tests of the system's behavior by separately seeding each 
individual node. Of these, 7865 contagion events were 
successful leading to long term non-zero infection levels. 
Figs. show three example system evolutions for 

the same network, differing only by seed node. In each 
case, the early dynamics of <fi have a chaotic appearance 
with a broad period three pattern, but then all collapse 
sharply to (apparently, see below) different kinds of peri- 
od 2 behavior. 

The time to collapse t c varies greatly, with Fig. [3p 
showing an example requiring over 30,000 time steps and 
over 3xl0 8 individual node updates. In Fig. [3f>, we pro- 
vide the complementary cumulative distribution for col- 
lapse times. The semi-log scale indicates that an expon- 
tial decay for i c <20,000 covers the majority of cases. A 
set of exceptional dynamics last for much longer, with 
one group ending at around t c ~ 30,000, and another 
isolated set distributed around t c ~ 25,000. Even after 
the dynamics collapse, we see that on the order of 10% of 
all nodes change state in each time step. We find similar 
behavior in for other values of fc avg and for larger systems 
of size N = 10 5 . 

We turn lastly to the behavior after the collapse. We 
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FIG. 3: Comprehensive survey of the possible Limited Imitation Contagion dynamics for a single fixed random network 
(fcav g =30) with its A?=10 4 nodes having fixed deterministic response functions. We ran N simulations, activating each node as 
a single seed and allowed the fully deterministic dynamics to run until collapse. (A-C): Examples of dynamics for three systems 
showing widely varying collapse time t c and different post-collapse dynamics. The inset in A shows one cycle of the underlying 
period 8 behavior; the post-collapse dynamics in B and C have periods 24 and 8. (D): Complementary cumulative distribution 
of collapse times. (E): Invariant densities; (F): Histogram of post-collapse microperiods; (G): Histogram of post-collapse 
macroperiods. 



call the period of <£> the system's 'macroperiod', and the 
periods of individual nodes 'microperiods'. Fig. 03 shows 
the invariant density \ for all seeds post collapse, ordered 
by x's variance (contagion events that fail are not includ- 
ed). The gray scale indicates frequency of the invariant 
density, and the blurring shows that higher order peri- 
ods underlie the apparent macroperiod 2 behavior. For 
example, the macroperiods of the time series in Figs.[3]f\.- 
C are 8, 24, and 8. We also note that the tent map's fixed 
point of 2/3 (horizontal line in Fig. OS) is not the center 
of the invariant densities, and that this is perhaps a finite 
size effect. 

In Fig. [3}? and G, we provide complete histograms of 
all post-collapse microperiods and macroperiods (again 
ignoring failed spreading events). We see that micrope- 
riods such as 1, 2, 4, 12, 24, and 48 form a dominant 
envelope, and while there are occasional instances of 
microperiods 3 and 5, many odd- numbered microperi- 
ods are absent. We find that the most common resultant 
macroperiods are 4, 8, 12, and 24; that a pure macrope- 
riod 2 is relatively rare (< 0.1%), and that the largest 
observed macroperiod is 120 (5 instances). Again, all of 
these outcomes are deterministic, depending only on the 
choice of initial seed. 

To conclude, in abstracting from a real world prob- 
lem, we have constructed a rich networked-based model 
of Limited Imitation Contagion that allows us to deeply 
explore the effects of increasing fixedness in a well-defined 



fashion. Some next steps would include, via simulation 
and analysis, exploring the long term behavior for the 
D-F class when there is no collapse 18j, investigating 



how collapse time distributions scale with TV for the D- 
F class; similarly testing other contagion mechanisms 
across systems of increasingly fixed microscopic struc- 
ture; and examining these Limited Imitation Contagion 
models on real networks and classes of naturally occur- 
ring ones such as scale- free networks [2fj| . 
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